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/\ ' We calculate the kaon S-parameter in quenched lattice QCD at /3 = 6.0 using Wil- 

5t ! son fermions at k = 0.154 and 0.155. We use two kinds of non-local ("smeared") sources 

for quark propagators to calculate the matrix elements between states of definite momen- 
tum. The use of smeared sources yields results with much smaller errors than obtained 
in previous calculations with Wilson fermions. By combining results for p = (0, 0, 0) and 
p = (0,0, 1), we show that one can carry out the non-perturbative subtraction necessary 
to remove the dominant lattice artifacts induced by the chiral symmetry breaking term 
in the Wilson action. Our final results are in good agreement with those obtained using 
staggered fermions. We also present results for S-parameters of the A/ = 3/2 part of the 
electromagnetic penguin operators, and preliminary results for Bk in the presence of two 
flavors of dynamical quarks. 
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1. Introduction 

Present calculations of weak matrix elements in the quenched approximation with 
Wilson fermions suffer from two main sources of error: (i) the signal is poor and (ii) 
there are large 0{a) corrections due to lack of chiral symmetry [0 0. In this paper 
we investigate the calculation of the matrix elements of four-fermion operators between 
pseudoscalar states, and in particular Bk- To improve the signal we calculate the 3-point 
function by sandwiching the operator between kaons produced by smeared sources. This 
trick has been used to obtain very accurate results with staggered fermions ^j. In order 
to reduce the 0{a) artifacts we use a momentum-subtraction technique similar to that 
tried earlier by the ELC collaboration [Q. We find that the combined method reduces the 
statistical errors for all four-fermion operators we have looked at, and allows us to perform 
non-perturbative subtractions for removing two of the three chiral symmetry violating 
terms in Bk- 

The 0{a) corrections arise due to mixing between operators of different tensor struc- 
ture induced by the explicit chiral symmetry breaking term introduced by Wilson to re- 
move lattice doublers. In principle this mixing can be calculated in perturbation theory, 
but there are large non-perturbative effects at values of g used in lattice calculations. 
There are two approaches to improving the situation: one is to work with an improved 
action so that the mixing occurs at 0{g'^a) and 0{a'^) rather than at 0{a) [Q, and the 
second is to devise non-perturbative methods to subtract off the lattice artifacts. It is 
likely that the eventual solution will be a combination of the two methods. To this end 
we demonstrate that the calculation of matrix elements within states of definite lattice 
momentum works for p = (0,0,0) and (0,0,1), and furthermore that one can reliably 
carry out a non-perturbative subtraction using these two values of momentum. We use 
the kaon S-parameter as the testing ground for two reasons: (a) there are very accurate 
results available using staggered fermions on the same set of lattices against which we may 
compare our results, and (b) there is no mixing with operators of lower dimension. 

To make our non-perturbative method work we need two kinds of hadron source: one 
that produces hadrons with zero momentum and the other that couples to all momenta. 
We construct zero momentum hadron correlators using wall source quark propagators, 
while the Wuppertal source propagators yield hadron correlators that have overlap 
with all momenta. We have shown in Ref. that these two kinds of correlators yield 
reliable signals for both the amplitude and the mass extracted from 2-point correlation 
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functions. That paper describes in detail the lattices used in the calculation and details of 

the quark propagators and hadron correlators. It also contains results for hadron masses 

and decay constants obtained from 2-point correlation functions. We use 35 lattices of size 

16"^ X 40 at /3 = Gg~^ = 6.0 with quark propagators calculated at k = 0.154 and 0.155. The 

two values of k correspond to kaons of mass Mk — 700 MeV and 560 MeV respectively, 

using a~^ = 1.9 GeV. 

The most accurate results for Bk at /3 = 6.0 have been obtained with staggered 

fermions [Q: 

„ ^r 0.70 ±0.02 : (16^ X 40 lattices) . . 

^ ~ \ 0.70 ± 0.01 : (243 X 40 lattices) ' ^ ' ' 

There are two previous estimates of Bk with Wilson fermions at (3 — 6.0. The results of 
Bernard and Soni are [0: 

„ ^ r 0.83 ±0.11 ±0.11 : (16^ X 40 lattices) , . 

^ ~ \ 0.66 ± 0.08 ± 0.04 : (24^ x 40 lattices) ' ^ ' ^ 

and that of the ELC collaboration is [jl| : 

Bk = 0.81 ± 0.16 ± 0.06 : (10^ x 20 x 40 lattices) , (1.3) 

where the second error is an estimate of the systematic error due to the subtraction of 
the bad chiral behavior. In all calculations the lattice kaon consisted of two almost de- 
generate quarks (the ratio ms/rriu < 3). The above results were obtained after interpo- 
lation/extrapolation of the lattice data to a kaon mass of 495 MeV. The large spread in 
these numbers and the systematic errors due to bad chiral behavior induced by the Wilson 
term underscore the need for further improvements and new methods. 

We also calculate the S-parameter for the A/ = 3/2 part of the electromagnetic 
penguin operators O7 and Og. Previous calculations with both Wilson M H, and staggered 



fermions [^ show that reliable results for the matrix elements of these LR operators can 
be obtained in lattice calculations and that the vacuum saturation approximation (VSA) 
provides a good estimate, i.e. Bj^ = 1.0 ± 0.1. Our estimates are 0.89(4) and 0.93(5) 
respectively, and we find that the dominant contribution to the matrix elements of both 
the LR operators and their VSA comes from the pseudoscalar (g) pseudoscalar (V) part of 
the 4-fermion operator. Our data show that matrix elements of V are larger by a factor of 
10 or more than other tensor structures and that the 2-color loop contraction is roughly 
three times larger than the 1-color loop. Furthermore, as the operator V is not suppressed 



in the chiral limit, we believe that VSA will be a good approximation in cases where the 
operator or its fierz transform contains V at tree level. 

This paper is organized as follows: in Section 2 we review the problem induced by 
the Wilson r term and our partial solution for subtracting lattice artifacts. In Section 
3 we describe the lattice methods and in Section 4 we present our results for Bk- We 
make a comparison with earlier results obtained with both Wilson and staggered fermions 
in Section 5. Section 6 presents preliminary results for Bk with two flavors of dynamical 
quarks. The analysis of the LR operators is given in Section 7 and we end with conclusions 
in Section 8. 



2. Bk and the problem of bad chiral behavior 

Weak interactions give rise to mixing between the K^ and K^. The relevant oper- 
ator in the low energy effective weak Hamiltonian is the /S.S = 2 four-fermion operator 
{s'-y ^Ld){s'-y ^Ld) , where we use the notation L = (1 — 75) and i? = (1 + 75). The value of 
the matrix element of this operator between a K^ and K^ at a typical hadronic scale is 
severely influenced by strong interaction effects. It has become customary to parameterize 
this matrix element by the kaon S-parameter, Bk, which measures the deviation from its 
value in the VSA: 

{K^\{s^M){-s^^Ld)\K'') = ^ImIBk, (2.1) 

where () indicates a trace over the spin and color indices. The normalization used for 
the decay constant is such that /^r = 132 MeV. If the VSA is exact then Bk = 1- To 
calculate Bk from flrst principles we must turn to non-perturbative methods such as the 
lattice. Our lattice calculation of Bk uses Wilson's formulation for fermions. The inherent 
violation of chiral symmetry in this approach leads to technical difficulties which we now 
review. 

To begin with, note that {s'-y^Ld){s'-y^Ld) is a special case of the operator 

0+ = ^((^i7^L^2)(^37^^V'4) + (2 ^ 4)), (2.2) 

with '01 = -03 = s and '02 = '04 = d. The significance of this is that with a chirally invariant 
regulator (9_|_ is multiplicatively renormalized. With Wilson fermions, however, this is not 
the case: there is mixing of this LL operator with other tensor structures in addition to 
an overall renormalization, and this complicates the definition of a lattice operator with 



the desired continuum behavior. In perturbation theory, the corrected operator has been 
calculated to 1-loop in Refs. |jTl[] and fl^: 



_|iz,(,„,))or+i|i 



0+°°' = (1 + tI;3^+('-. «''))0+" + jk^ r''Z-(r)(Ol^'' + Ol-^ + Of) (2.3) 

where 



OV = - "^'32^"' [(V - ^) + (2 - 4)] , (2.4) 

^r=-Y^[('5-^)+(2-4)], 

and A^ = 3 is the number of colors. We have used a condensed notation for the allowed 
Lorentz tensor structures: 

S = (-01^2) (^3V'4), 

V= (V'i7mV'2)(^37mV'4)- 

^= 5Z(^l^/^^^2)(V'3(7Mt'^4), (2.5) 

A = (^i7m75V'2) (^37^75^4), 

V = (^175V'2)(V'375'04), 

where 7/^,75 are hermitian and a^j^ = (7^71/ — 7i/7^)/2. We note that the Fierz transform 
eigenstates appearing in Eq. ( pT3|) are only (V+^), |(V— ^)±(5— P) and (iS+T+P); there 
is also no mixing between the fifth eigenstate of the Fierz transformation {S — ^T + V) and 
the operator Oj^ at 1-loop. There is no mixing with lower dimensional operators, for the 
simple reason that there are no AS = 2 operators of lower dimension. We shall henceforth 
denote the perturbatively corrected {s'y^Ld){s'yfj,Ld) operator (cf. Eq. (|2.3| )) by O. 

The renormalization coefficients for Wilson parameter r = 1 are given in Table 1 in 



three schemes: the dimensional reduction (DRED) used by Altarelli et al. [0 and Mar- 



tinelli |jTT[], as well as the "naive" dimensional regularization (NDR) and the dimensional 



reduction {DR{EZ)) scheme used by Bernard et al. in |jT^. A detailed description of 
DRED and NDR schemes and their relative advantages and disadvantages is given in 



Ref. [|I^ . We tabulate the relevant results in order to provide easy reference, and to allow 
the reader to make a rough estimate of the magnitude of the scheme dependence. All 
our results are given in the DRED scheme, except when we compare raw lattice numbers 



against those in Ref. [jT5|, in which case we use DR{EZ). 



For each of the four-fermion operators, S, V, T, A, and V, there are two distinct 
contractions with the external states. In the first each bihnear is contracted with an 
incoming or outgoing kaon corresponding to two spin and two color traces. We label these 
contractions by V^, iS^, V^, A"^ and T^. The other contraction consists of a single spin and 
color trace which we Fierz transform to two spinor loops. We label them by V^, S^, V^, 
A^ and T^, since they have a single color trace. We will find it useful to further split the 
V, A and T terms into their space and time components, and denote these components 
by subscripts s and t respectively. This notation is similar to that used with staggered 
fermions and will facilitate later comparison of results for individual operators between 
the two formulations. 

In order to extract Bk, we calculate, at non-zero momentum transfer, the matrix 
element 

Mk{p) = (K^ip)\dip)\K\p = 0)). (2.6) 

In chiral perturbation theory A4k{p) behaves as ~ ■JkPk ' Pk^ where 7^^ = 8/3 /^S^, 
and pk and pp^ are the on-shell four-momenta of the external states, so that pk ■ PR = 
M K yM^^+ip)^ . Unfortunately, on the lattice with Wilson fermions chiral symmetry is 
explicitly broken and the expansion becomes 

Mk{p) = a + (3Ml + (7 + 7K W ■Pk + ---- (2.7) 

Here the terms proportional to a, /3 and 7 are unphysical contributions arising due to the 
r-term in the Wilson action, and suppressed by one power of the lattice spacing a. Similar 
formulae hold for each individual spin-color term described above, and apply to both the 
on-shell {{K^\d\K^)) and the off-sheU {{{)\d\K^ K^)) matrix elements. 

Using O should reduce the lattice artifacts, but it will not eliminate them completely 
because it is only an approximation to the operator with the desired continuum behavior. 
In particular, the coefficients Z contain terms of 0{g^) and higher that have not been 
calculated, and, more importantly, as previous calculations have shown, there are large 
non-perturbative effects. We therefore require non-perturbative methods to isolate the 
physical coefficient 7^. 

The most troublesome of the lattice artifacts is a. Failure to correctly subtract this 
contamination will mean that Bk will diverge in the chiral limit. To eliminate this, it is not 



necessary to work at non-zero momentum transfer. One can simply calculate M.k{p = 0) 
at different values of n (that is, different values of Mk) and take a difference, leaving 

MK{^^,p = 0) - A^k(/^,P = 0) = (/3 + 7 + 1k){MI - Ml). (2.8) 

To remove /3 one takes the difference of the on-shell and off-shell matrix elements. This 
method has been used in Refs. |^] and 0, and suffers from the lack of control over final 
state interactions between the kaons in the off-shell amplitude. A review of the status of 



previous results is given in Refs. [jl^] and []17 



We advocate using momentum subtraction which eliminates both a and /? at a fixed 
value of Mk, using only on-shell quantities. For example, by calculating the matrix element 
of O for two different values of p and taking the difference one gets 

Mk{p) - Mk{Q) = (7 + 1k)Mk{E{p) -Mk) + .... (2.9) 

In practice we calculate 

E(p) B^(p) - Mk Bj^(0) , , , , 

at each value of k, where by Bj^{p) we mean the ratio of the matrix element to its VSA 
value, both calculated on the lattice at appropriate momentum transfer. 

This momentum subtraction scheme does not eliminate the lattice artifact 7. Fur- 
thermore, in working at non-zero momentum there is a danger that higher order terms 
(for example quartic in momenta) omitted in Eq. (^) may become significant. Therefore 



in order to compare the lattice result with experiment we have to make the following as- 
sumptions: (i) using the perturbatively improved operator O makes 7 negligible; (ii) the 
terms of order p^ and higher that we have neglected in the chiral expansion do not have 
large coefficients. The first assumption is expected to become more reliable on using an 
improved Wilson action, while the second will come under better control as calculations are 
done on larger spatial lattices since then the gap between lattice momenta will decrease. 
At this stage the only justification for these assumptions is the a posteriori agreement of 
results with those obtained using staggered fermions. 
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3. Methodology 

Our method for calculating B^ requires that we double the 16^ x 40 lattices in the time 
direction, so that they are of size 16^ x 80. On these doubled lattices we construct hadron 
correlators such that the correlator on time slices 1-39 is the forward moving particle 
with the source at time slice 0, while the correlator on time slices 79-41 is the backward 
moving particle with the periodically reflected source on time slice 80. To calculate matrix 
elements we insert the operator between these "forward" and "backward" moving particles 
on the original 16'^ x 40 lattices. 

In practice, we divide the correlators for the various matrix elements by the product 
of kaon correlators, so that we directly obtain the S-parameters for the various operators 

O 

^ ^3 _{K^{p)\O\K^{p = 0)) 
"" 8{K0{p)\A,\0){0\A,\K0{p = 0)y 

We can select the kaon momenta by our choice of source and by inserting momentum into 
the operator O. The statistical errors are reduced because we can average the operator 
location over a time slice of the lattice. Away from the sources, only the lightest state 
contributes to the correlators, and we should flnd a time-independent plateau giving Bo- 
This method is similar to the one we have used successfully with staggered fermions [Q. 

The physical picture of the process for calculating matrix elements using smeared 
sources is as follows: a wall source at t = produces zero momentum K^ which propagates 
for a time t, at which point the operator inserts momentum p, and the resulting K^ with 
momentum p then propagates the remaining (Nf-t) time slices until it is destroyed by a 
Wuppertal source. Three factors are essential for our method to work: 
(i) The wall source creates only zero-momentum kaons; otherwise there is contamination 

from matrix elements of kaons with other momenta. 
(ii) The Wuppertal source has signiflcant overlap with the lowest few momenta allowed 

on the lattice, 
(iii) For matrix elements involving p ^ kaons, we must ensure that there exists an 

overlap region for the kaons where a plateau can be observed in the S-parameter 

signal. Thus it is essential that the signal for the zero-momentum kaon produced by 

the wall source extends across the lattice to the region where there is a signal for the 

non-zero momentum kaon produced by the Wuppertal source. 



In Ref. 0] we showed that these conditions are satisfied by the Wuppertal and waU cor- 
relators, when we use p = (0,0,0) and p = (0,0, 1). Furthermore, there are a number of 
consistency checks we make: 

(1) The p = matrix element is calculated three difi'erent ways; (a) using wall sources on 
both sides, (b) using wall source on one side and Wuppertal source on the other, and 
(c) using Wuppertal sources on both sides (in this case there is a small contamination 
from the p = (0, 0, 1) state). 

(2) We use two kaon source operators: 75 and A4. The plateau in each individual Bq is 
reached from opposite directions for these two. The two results should converge to 
the same value. 

As shown in Tables 2a-3b, these checks are satisfied by our data within the statistical 
accuracy. We also find that the S-parameter for the operators J^ and V^ are within a 
factor of two of their VSA values. As in the case of staggered fermions, the final value of 
Bk is obtained after a large cancellation between the A and V components showing that 
VSA is not a good approximation. 

4. Results for Bk 

Our final lattice result at a given value of k and p is obtained from the perturbatively 
improved combination (using the convention that all four quarks have distinct fiavor labels 
so that each term has just one Wick contraction) 

Zisi(p-)=('l + ^Z+(r,a^) 



g2 ^2^*^^^ 



167r2 12 



V^ + V^ + ^^ + A^ 
(265^ + 25^) - (18P^ - 6P2) (4.1) 

+ 4(Ti + T^) + (V^ - A^) - 11(V2 - A^) 



For simplicity, we have here used the operator symbol to denote its S-parameter. The 1- 

loop perturbative results for the renormalization constants Za and Z_|_ are given in Table 

2 
1. Note that the finite part of the renormalization factor (1 + ^5:^-^+) is largely canceled 

by Z\ for g^<.l. 

We give the results for the individual Bo (without any g"^ corrections) at k = 0.154 

and 0.155 in Tables 2a-3b. In all cases we find that the signal in the ratio of correlators 

is significantly better with the operator 75 as the kaon source than with the operator ^4, 



even though the two sets give consistent results. As an example, we show a comparison 
of the two signals in Figs, la and lb. Our final results therefore use the 75 numbers. 
We point out that in case of non-zero momentum transfer, the signal for Bq only exists 
closer to the Wuppertal source (at time slice "40") than the wall source (at time-slice "0"). 
This is because the signal in the p = (0, 0, 1) kaon correlators only extends for about 20 
time-slices. The overall quality of the signal for B^ (with g'^ = 1, fia = 1.0) is shown in 
Figs. 1, 2, 3 and 4 at various values of k and p. 

Tables 2a-3b show that all individual Bq, except for -B^^, increase by a factor of 
about 2 between k = 0.154 and 0.155. This increase is largely due to the change in VSA, 
i.e. the factor f'^M]^ decreases by approximately 1.9 between the two k values 0. This 
shows that at these values of k, the lattice matrix elements are dominated by the constant 
term a. 

The contribution of the mixing terms to B^ can be large only if the matrix elements 
are large, since the perturbative mixing coefficient is ~ 0.005 for g^ = 1. The data show 
that the largest matrix elements are of the operator V; however their net contribution to 
B^ is very small, since V^ ~ 3V^ (approximate VSA). Both T^ and V^ are close to zero. 
The next largest contribution comes from 4T^, which is partially canceled by 265^ + 2S^. 
The net result of these features in the data is that the contribution of mixing terms to 
Bx is in fact of the order of a few percent. Unfortunately, since the unphysical term 7 in 
Eq. ( |2.7|) also gets contributions from the diagonal operators, the small value of the mixing 
terms does not provide a bound on 7. 

Given B^{p) we calculate Bk and the errors using Eq.( p.lO| ) two ways: (1) for each 
jackknife sample we first perform the momentum subtraction and then the mean value and 
the error are obtained as the jackknife estimate over the 35 samples, and (2) we construct 
the four quantities needed in Eq. (|2.10| ) independently along with their errors, and obtain 
the final error estimate assuming that the individual estimates are uncorrelated. Our 
quoted results use the first method, but we note that both the methods yield consistent 
estimates. 

We have calculated Bj^{p = (0, 0, 0)) three different ways: using Wuppertal- Wuppertal 
{S — S), Wuppertal- wall {S — W), and wall- wall {W — W) correlators. For example, our 

9 



results, using g'^ = 1.0 and ixa = 1.0, are: 



0.38(7) 

B^{k = 0.154) = { 0.37(4) 

0.37(4) 



S-S 
S-W , 

w-w 

(4.2) 
0.10(11) -.S-S 

B^{k = 0.155) = { 0.11(7) : S-W . 

0.13(7) : W-W 

The consistency of the data suggests that the contamination in the S — S result from 
higher momentum kaon states is at most a few percent. Since only the S — W correlators 
give results at both p = (0, 0, 0) and (0, 0, 1), we shall henceforth quote results for Bj^{p) 
obtained from these. 

In order to extract the continuum result for Bk we must choose the values of both g^ 



and jja to use in Eq. (|4.1|) . Lepage and Mackenzie have suggested [jT^ that perturbation 
theory is much better behaved if one uses the coupling constant in a continuum scheme 
such as MS, instead of the bare lattice g'^. They also give a prescription for choosing 
the appropriate scale of the coupling constant. In general this scale will differ for the 
various operators that mix with (9_(_ in Eq. ( [4.1|) . To simplify the calculation we take all 
the scales, and thus all the coupling constants, to be the same, i.e. all of 0{Tc/a). Then 
the Lepage-Mackenzie prescription amounts to a replacement of the bare lattice g^ with 
an effective coupling g'^^ ^ ^■75g'^ ai (3 = 6.0. To study the dependence on g'^^f, we use 
four different values, g'^g = 0.0, 1.0, 1.338 and 1.75. It is important to realize that only a 
2-loop calculation of the perturbative coefficients can test whether a given choice for g'^^^ 
is reasonable. Such calculations that have been done to date support Lepage-Mackenzie 
prescription for choosing g'^^ [ [T^ . 



The choice of fia is of a different character to that of g^. In physical matrix elements 
(e.g. that related to CP- violation in K^ — K^ mixing) B^ always appears multiplied by a 
coefficient function, such that the combination is independent of //. At leading order, the 
scale independent combination is 

SK = a.(^)-'/(''-'^^)i?K(^), (4.3) 

where Nf is the number of active flavors. In fact, Bk does have some dependence on 
/i, coming from the following sources. First, since we are using only the leading order 
expression for Z{fxa), Bk does depend on fx at non-leading order: dlnBK/dlnfx oc (/^(/u). 
This is likely to be a small effect, and it can probably be pushed to next order given the fact 
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that the two-loop anomalous dimension and one-loop matching coefficients are known. We 
say "probably" because it is possible that there are some residual subtleties with Wilson 
fermions associated with the mixing of 0+ with opposite chirality operators. A related 
source of fj, dependence occurs when fia differs greatly from unity: then higher order terms, 
proportional to [g'^ \n{fxa)]'^ , which are not included in Eq. ([4.1|), become large. What is 
happening is that the leading logarithms, which have been summed into the coefficient 
function, are partially incorporated into the perturbative coefficients. Once again, one can 
probably take these into account knowing the anomalous dimension to 2-loops, or finesse 
the problem by taking |Ua ~ 1. Finally, we are calculating the lattice result in the quenched 
approximation, for which the number of active flavors is zero, while we wish to match to 
the full theory with Nf active flavors. This introduces a small /i dependence. 

Our emphasis in this paper is on improving methods for calculating Bk, and not on 
extracting final numbers for B^- Thus we choose to quote our results for a variety of 
values of fia so as to allow others some flexibility if they wish to use our numbers. We use 
/ua = 1.0, TT, and I.Ttt. We have a slight preference for /xa = tt, since then the continuum 
and lattice cut-offs are matched. 

Table 4 shows the sensitivity of the results to the choices of parameters, for both 
B^ and the momentum subtracted Bx- There is a significant variation of the results 
with g'^^ and fxa. For a fixed value of /ua, the lattice results for both Bjj^{j> = (0,0,0)) 
and B^{p = (0,0,1)) increase as a function of (7^^ due to the increased contribution of 
the mixing operators. For fixed g'^g, an increase in fxa decreases the contribution of the 
diagonal operators (note that B^{p = (0,0,0)) at k = 0.155 is insensitive to changes in 
fia because the diagonal contribution happens to be almost zero there). 

As for Bk after momentum subtraction, the estimate decreases by 10 — 20%, at both 
values of quark mass, between g'^^^ = 1.0 and g'^^^ = 1.75, for fixed fxa. It turns out that 
almost all the variation comes from the diagonal renormalization constants (for example, 
the ratio of (1 + {g'^/16n^)Z+) to Z\ changes from 0.93 to 0.8), and not from operator 
mixing. As our present best estimates for Bk we quote the results at g'^^ = 1.75, and use 

fxa = n: 

BKiti = 0.154) = 0.68(22) , 

(4.4) 
Bk{k = 0.155) = 0.57(23) . 
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5. Comparison with previous calculations 

The staggered fermion results for Bk are statistically the most accurate and have 
the correct chiral behavior. At /3 = 6.0, the kaon mass roughly matches for staggered niq = 



(0.02 + 0.03) and k = 0.154, and for staggered ruq = (0.01 + 0.02) and k = 0.155 0. Thus 
we can compare the corresponding data for B^- This is a particularly good comparison 
because the two calculations have been done using the same set of background gauge 
configurations. Using gf^ = in the 4-fermion renormalization constants, the staggered 
results are 0.76(1) and 0.72(2) to be compared with our numbers 0.83(21) and 0.77(19) 
respectively. A striking feature is that the errors with Wilson fermions are a factor of 10 or 
more larger. The data in Table 4 indicate that the errors in individual i?^(p) are larger by 
about a factor of four, and the remaining factor comes from the momentum subtraction. 
The staggered results were obtained by making fits without including the full covariance 
matrix, and if we do the same for Wilson fermions then the errors in individual B^{p) are 
reduced by a factor of about two. Thus, at the level of the signal in the correlators, the 
smeared sources work almost as well for Wilson fermions as for staggered fermions. It is 
the process of momentum subtraction that leads to a significant increase in the error. 

It could be argued that one should actually compare staggered numbers with our pre- 
ferred results using f/gg- = 1.75 and na = tx, i.e. 0.68(22) and 0.57(23) respectively. The 
rationale for this is that a good choice of g^^^ and fia aims to reduce the effects of artifacts 
a, /3 and 7. Recall that a and jS are absent for staggered fermions while a 7 like term exists 
due to operator normalization. With this choice the Wilson fermion estimates lie system- 
atically below the staggered values. The perturbative corrections for staggered fermions, 
though small, have not been included in the published results of B^ [^ • Including them 
would decrease Bk by ~ {2g'^g)%. The remaining difference could be due to the difference 
in the 0{a) corrections for the two fermion formulations, e.g. the artifact 7. In any case, it 
is not clear whether the difference is significant, given the size of the errors in the Wilson 
results and the large variation with g'^^ and fxa. 

We can also make a direct comparison with results obtained by Bernard and Soni using 
Wilson fermions [0 . They have calculated Bj^ {p = 0) using the perturbatively improved 
operator on a subset of the same lattices (they used only 19 lattices as they skipped every 
other one in each of the two streams), and at the same two values of k. Their method of 
extraction of B^ is described in Ref. [|2|, and is different than the method we have used. 
Using (/2^ = 1.338 and fia = 1.7ti then resuhs are B^{p = 0) = 0.36(22) and 0.24(45) 
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at K = 0.154 and 0.155 respectively [^, to be compared with B^(p = 0) = 0.38(5) and 
0.11(7) obtained by us. Even after allowing for a factor of ~ \/2 due to statistics, it is 
clear on the basis of this comparison that our use of smeared sources has reduced the errors 
considerably. 

Though our results for Bj^ are more accurate compared to those obtained by other 
groups using Wilson fermions, our final results for Bx have larger errors (cf. Eq. ( |1 . 2| ) , 
( p..3| ) ) . The gain due to the use of smeared sources is compensated by the increase in error 
due to momentum subtraction. The advantage of using momentum subtraction is that it 
unambiguously removes lattice artifacts a and /3. Also, numerical errors in Bj^{p = 1), 
as well as the contribution of quartic terms in the chiral expansion, should decrease when 
using a larger lattice due to the decrease in the value of lattice momenta. 

One further qualitative comparison that we can make is for the S-parameters (without 
perturbative improvements but after momentum subtraction) of the individual space/time 
and l-loop/2-loop components of the four-fermion vector and axial operators, with the 
corresponding results obtained using staggered fermions . Such a comparison is possible 
because, as discussed above, the effects of operator mixing are small. This comparison 
provides information on the reliability of the momentum subtraction procedure for Wilson 



fermions. Furthermore, as explained in Ref. |21|], the chiral behavior of Sy and i?_4 is 
known; both are expected to increase in magnitude with decreasing quark mass due to 
chiral logarithms and finite volume dependence, and can therefore provide a sensitive test 
at small quark masses. The results of our comparison are shown in Table 5. Though the 
errors in the results with Wilson fermions are much larger, it is reassuring to see that the 
central values are in good agreement. In fact the agreement is far more impressive than 
the errors would naively lead us to believe. We need to perform calculations at more values 
of K and P to confirm this favorable behavior. 

6. Results with two flavors of dynamical fermions 

We have also estimated Bk, using the same methodology as above, on 16^ lattice 
configurations generated with two fiavors of dynamical Wilson quarks. The details of 



these lattices are given in Ref. |2^. The kaon mass at /? = 5.5, k = 0.159 and 0.160 is 
roughly Mk = 860 MeV and 650 MeV respectively; and at /3 = 5.6, k = 0.156 and 0.157 it 
is about 1050 MeV and 820 MeV respectively. This calculation of Bx has been done only 
with kaons created with a 75 operator. We find that a time interval of 16 is large enough 
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to give a stable plateau over about 6 central time-slices for both p = (0, 0, 0) and (0, 0, 1) 
correlators, as illustrated in Figs. 5 and 6. The behavior of individual Bq terms is very 
similar to the quenched numbers shown in Tables 2a-3b. 

Before presenting the final results, we first discuss a technical drawback of this cal- 
culation. On these lattices only Wuppertal source correlators are available, so with our 
method some contamination from higher momentum kaon states is present in the data. For 
momentum transfer p = (0,0,0), the largest contamination comes from the propagation 
oi p = (0, 0, 1) kaons across the lattice. (The contamination in the p = (0, 0, 0) -^ (0, 0, 1) 
data comes from the presence oi p = (0, 0, 1) -^ (0, 0, 2) terms.) This contribution is 
suppressed by two factors: the exponential suppression due to the extra energy of the 
p = (0, 0, 1) state, and the square of the ratio of amplitudes for creating a, p = (0, 0, 1) 
kaon versus a p = (0, 0, 0) kaon by the Wuppertal source. These factors increase as the kaon 
mass decreases. They are similar on our four sets of lattices. At /3 = 5.6 and k = 0.157, 
our estimates are ?« 10 and (1.5)^. There is also an enhancement factor because the matrix 
element between higher momentum kaon states is larger. We estimate this factor using 
VSA to he {E{p = 1)/Mk)'^ ~ (1.5)'^. These three factors combine to increase the result 
for Bk by roughly 10%. We note that in the case of the quenched lattices, having 40 
time-slices reduced this contamination to the level of a few percent. This is evident on 
comparing the S — S and the S — W or the W — W results in Tables 2a-3b. 

We again calculate B^ for three values of the effective coupling: 0.0, g^ and 1.75g^. 
The lattice scale on individual lattices (without extrapolation to the chiral limit) is not 
well defined, and we simply set fia = 1. The results are shown in Table 6. As explained 
above, the results for Bk are likely to be ~ 10% larger due to contamination from higher 
momentum kaon states. In addition, Bk has to be extrapolated to the physical kaon mass. 
Thus, the only conclusion we can draw is that the quenched and dynamical results are in 
qualitative agreement for quarks masses in the range nis < ruq < 3ms . 

7. S-parameter for the Left-Right electromagnetic penguins 

There are two additional 4-fermion operators that we analyze using the data in Tables 
2a-3b. These are the A/ = 3/2 part of the left-right electromagnetic penguin operators O7 
and Os- They alone contribute to the imaginary part of the 1 = 2 amplitude and therefore 
give the dominant electromagnetic contribution to e'/e. A knowledge of their S-parameters 



is phenomenologically important as discussed in Ref. [^. Taking just the A/ = 3/2 part 
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of the operators simplifies the numerical calculation as the "eye" contractions cancel in 
the flavor SU(2) limit. 

In principle one would like to calculate the matrix elements of the penguin operators 



O7 = (sai^Lda) [(ub'j^.Rub) - -(db'j^.Rdb) - -{sb'^^^Rsb)]., 



Os = {Sal^iLdb) [{Ub'^^^RUa) - -{db^^Rda) - -{Sb'J^.RSa 



(7.1) 



(7.2) 



between a K~^ and a 7r+. Instead, we calculate the A/ = 3/2 part given by the operators 

Cy^^ = (Sa'J^^Lda) [{Ub'J^^RUb) " (db'J^.Rdb) ] + {Sa'J^^LUa)iUbl^Rdb), (7.3) 

^8^^ = (sa'j^j.Ldb) [{ub'y^.Rua) - {db'jf.Rda) ] + {sa'y^^Lub){ub'y^^Rda)■ (7.4) 

Note that the overall normalization is unimportant as it cancels in the S-parameters. The 
1-loop perturbatively corrected versions of these operators have been calculated in Refs. 
plj] ||12|| , and are linear combinations of the operators labeled Oi and O2 therein. The 



M7{p)=[l + j^Ziir,a^i) 



+ 



g" {Z2-Z^ 



matrix elements of these corrected operators between a K^ and a 7r+ are, in the flavor 
SU(2) limit, 

o^ \ r 

2p2 _ 252 + V^ - ^^ 
4pi _ i2p2 ^ ;l65^ - 165^ 

- 7V^ - iiv^ - 9^^ - hj? - er^ + 2r2 

and 



167r2 3 

g2 ^2^*^^^ 



1671^ 



12 



(7.5) 



MM=\\^^Z2{r,a^x) 



gi r^z*{r) 



167r^ 



12 



2p2 _ 252 ^ yl _ ^1 

18P^ - 38V'^ + US^ - 265^ 

- 5V^ -v^ + A^ - 3A^ - ler^ 



(7.6) 



where, if necessary, we have made a spin Fierz transformation to recast all the terms as 
two-spinor loops. The corresponding VSA contractions are 



,VSA 



Mr^'ip) 



'-ZlV'^ -ZlA^ 



(7.7) 
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and 

iVSA, 



M^^ip) 



2ZlV^ - -ZlA'^ 



(7.8) 



The S-parameters are the ratios of, for example, the matrix element of O7 to its VSA. 
We evaluate these in the SU(3) limit, i.e. degenerate w, d and s quarks. The 1-loop values 
of the renormalization constants for these LR operators are also given in Table 1. Note 

3 /2 

that in the DRED scheme the operator O^ does not mix with the scheme dependent 
operator O of Ref. [|T^. It is for this reason that we choose this scheme, although, our 
analysis shows that the results are only weakly scheme dependent. 

In the chiral limit these matrix elements are expected to behave as c + dm^ + . . .. 
There are 0{a) corrections in the coefficients c and d due to the lattice discretization. At 
present the only way to reduce these is to use an improved action and/or work at weaker 
coupling. In this study we do not have any control over these corrections and we simply 
give the lattice results for the Wilson action. 

The quality of the signal is shown in Figs. 7 and 8. In the analysis of these LR 
operators we find that using the full covariance matrix produced estimates that are about 
Icr lower than the fit values shown unless we significantly decrease the range of the fit. The 
error estimates with and without using the full covariance matrix in the fits are essentially 
the same. This indicates that the data at different time slices is highly correlated and larger 
statistics is needed to reliably include the correlations. We choose to use the full range of 
the plateau in the fit and quote results obtained without including the correlations. 

As in the case of the LL operator, in order to quote a value for the S-parameters we 
have to specify the value of g^^ and iia used in the perturbatively improved operators. 
In Table 7 we quote results for a number of choices in order to give an estimate of the 
sensitivity of the results to variation in these parameters. The data show that this could 
be a 10% effect, so it is important to make a good choice of g^^. 

The data also show a small increase in the S-parameters as the quark mass is de- 
creased. Linearly extrapolating the g^^ = 1.75, fxa = n results to the physical kaon mass, 

our best estimates are 

Sj/^ = 0.89(4), 

3/2 (7-9) 

Sg/^ = 0.93(5). 

These values are slightly smaller than the numbers used by Lusignoli et al. [p3| in their 
analysis of e'/e; they used S7 = Bg = 1.0 ± 0.1. To make a complete determination of 
e'/e we need to calculate many other matrix elements, for example of the strong penguin 
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operators O5 and Oq, for which the lattice technology is still unreliable. For this reason 
we do not consider it opportune to repeat the phenomenological analysis. 

3 /2 

We can make a direct comparison of lattice results for (9-, with those obtained 
by Bernard and Soni on the subset of lattices described in the analysis of Bk- Their 
result, obtained using g'^^ = 1.338 and fxa = I.Ttt in the DR{EZ) scheme, is 0.965(41) at 
K = 0.155, to be compared with 0.971(50) obtained by us. This comparison suggests that 
for the matrix elements of LR operators, our method of sandwiching the operator between 
smeared sources is no better than using propagators from a single source point. On the 
other hand the fact that smeared sources yield a plateau over a large range of time-slices 
gives reassurance that one potential source of systematic error is under control. 



8. Conclusions 

We show that the calculation of the kaon S-parameters with Wilson fermions is sig- 
nificantly improved by the use of non-local quark sources. By using a combination of 
Wuppertal and wall source correlators, we demonstrate that the on-shell matrix elements 
can be calculated at non-zero momenta. 

By combining results aX p = (0,0,0) and (0,0, 1), we carry out a non-perturbative 
subtraction of the lattice artifacts in the calculation of B^- Even though we cannot 
take into account the artifact 7, our results are in good agreement with those obtained 
with staggered fermions. On the basis of this exploratory study we feel confident that 
the momentum subtraction procedure indeed works. To make further improvements and 
reduce the 0{a) artifacts one needs to repeat the calculations with an improved lattice 
action and on a larger physical lattice with smaller Pmin- 

We find a clean plateau in the data for the S-parameters of the LR electromagnetic 
penguin operators. The results show that VSA works much better for these operators. 
All the B-parameters vary significantly with the choice of g'^^ used in the perturbative 
renormalization coefficients. Our final estimates are given using the value advocated by 
Lepage and Mackenzie in Ref. [jl^, i.e. g'^^ = 1.75. 

The method of using the combination of Wuppertal and wall correlators can be ex- 
tended to study other 3-point correlation functions, in particular structure functions and 



form factors. This work is in progress. 
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K = 0.154 and for momentuni transfer p = (0, 0, 0). The data are obtained using operator 
75 as the kaon source. 
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Fig. lb. The same as in Fig. la, but using operator A4 as the kaon source. 
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Fig. 2. The same as in Fig. la, but for n — 0.154 and p = (0, 0, 1). 
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Fig. 3. The same as in Fig. la, but for k — 0.155 and p— (0, 0, 0). 
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Fig. 4. The same as in Fig. la, but for «; = 0.155 and p — (0, 0, 1). 
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Fig. 5. The ratio of correlators for the lattice parameter i?^((7gg = 1.071, fxa = 1) mea- 
sured on 16^ lattices generated with two degenerate flavors of dynamical Wilson fermions. 
The lattice parameters are f3 = 5.6, k = 0.157 and momentum transfer p = (0,0,0), and 
we use S75(i as the kaon source. 
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Fig. 6 The same as in Fig. 5, but for p = (0, 0, 1). 
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Fig. 7 S-parameter for the LR electromagnetic penguin operator Oj . Tlie data are 
obtained using g'^^ = 1.75 and fxa = tt. 
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Fig. 8 S-parameter for the LR electromagnetic penguin operator O^ . Tlie data are 
obtained using g'^^ = 1.75 and fxa = tt. 
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DR{EZ) 


DRED 


NDR 


Za 


1 - 15.796CfA 


O.bCpX 





Zp 


1 + CFA(61og(H - 21.596) 





-CpX 


z+ 


-50.174 -41og(H 


7 
3 


14 
3 


z_ 


-45.308 + 81og(H 


2 
3 


4 
3 


Zx 


-49.364 - 21og(H 


11 

6 


11 
3 


Z2 


-42.064 + 161ogM 


8 
3 


16 
3 


z* 


9.6431 









Table 1. Summary of 1-loop perturbative results for the various Z factors needed in our 
calculations in three different continuum regularization schemes. The two constants are 
Cp = 4/3 and A = g'^/{16n'^). The results in DRED and NDR schemes are the sum 



of those in DR{EZ) and the entries in their respective columns. These expressions are 
extracted from Refs. |TI|] |T^ and 0. The numerical results are taken from Ref. |l'4j . 
In the text all results are given in the DRED scheme used in Ref. |]Tl|], except when we 



compare raw lattice numbers against those in Ref. |15 . 
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Table 2a. The one color loop contribution to the lattice S— parameters for individual 
operators at k = 0.154. Each box shows the x^ for a correlated fit, the temporal range 
of the fit and the fitted value. The space and time components of the operators have 
been shown separately. The appropriate ratios of correlators have been calculated using 
two different kaon operators, 75 and A4, and using Wuppertal and wall sources. The 
notation, for example, is: PS'[1]VFP stands for the four-fermion operator with one unit 
of lattice momentum sandwiched between Wuppertal and wall source kaons each created 
with operator 75. All errors are calculated using the single elimination jackknife method. 
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1.9 


2.1 


2.4 


2.1 


1.9 


1.2 




8-32 


11-26 


10-30 


25-34 


10-30 


14-26 


10-30 


25-33 


Al 


-0.54(5) 


-0.52(4) 


-0.53(4) 


-0.31(3) 


-0.50(7) 


-0.49(4) 


-0.49(5) 


-0.25(4) 




3.4 


4.2 


1.3 


0.8 


2.8 


3.3 


0.7 


0.8 




13-24 


10-25 


15-25 


25-33 


12-24 


14-26 


16-24 


24-32 


Ai 


0.64(4) 


0.62(3) 


0.58(2) 


0.64(8) 


0.56(5) 


0.58(7) 


0.54(3) 


0.60(7) 




5.3 


3.6 


3.2 


2.0 


3.6 


3.9 


3.0 


1.5 




10-30 


12-28 


10-30 


28-34 


12-28 


14-26 


10-30 


23-35 


Tl 


0.05(2) 


0.06(2) 


0.05(1) 


0.03(1) 


0.05(2) 


0.05(1) 


0.04(2) 


0.04(1) 




7.1 


1.5 


2.3 


1.0 


3.7 


2.1 


4.1 


1.1 




14-26 


14-26 


10-30 


25-33 


8-31 


14-26 


6-34 


26-33 


%' 


0.07(2) 


0.05(1) 


0.05(1) 


0.03(1) 


0.03(5) 


0.04(3) 


0.06(2) 


0.05(1) 



Table 2b. The same as Table 2a, but for the two color loop contribution. 
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PS[0]SP 


PW[0]WP 


PS[0]WP 


PS[l]WP 


AS[0]SA 


Aiy[o]iyA 


AS[0]WA 


AS[1]WA 




2.6 


0.6 


0.5 


2.1 


1.9 


1.8 


0.9 


0.8 




10-30 


16-27 


10-30 


23-33 


12-28 


14-28 


12-30 


26-33 


v 


-7.08(97) 


-6.48(41) 


-6.49(39) 


-3.42(57) 


-5.98(71) 


-5.99(44) 


-5.58(42) 


-3.24(42) 




3.2 


4.4 


1.8 


1.4 


4.5 


3.4 


1.5 


2.3 




10-30 


14-27 


11-29 


28-33 


10-30 


14-28 


12-30 


26-33 


s' 


-1.01(18) 


-0.92(19) 


-0.94(10) 


-0.48(13) 


-0.85(14) 


-0.89(19) 


-0.85(10) 


-0.48(11) 




5.9 


8.7 


3.1 


1.2 


4.6 


10.6 


3.3 


1.1 




13-29 


13-27 


13-27 


25-33 


14-32 


10-30 


12-30 


27-34 


vi 


0.83(65) 


0.87(41) 


0.91(17) 


0.32(18) 


0.83(18) 


0.83(31) 


0.79(18) 


0.36(27) 




6.0 


2.9 


1.5 


1.1 


5.3 


1.5 


1.0 


0.5 




15-26 


11-29 


14-26 


26-33 


11-30 


14-28 


13-29 


26-33 


v^ 


0.38(15) 


0.45(9) 


0.42(5) 


0.22(6) 


0.42(9) 


0.41(8) 


0.40(5) 


0.23(7) 




4.4 


2.6 


1.6 


0.3 


0.7 


3.9 


2.2 


0.8 




10-30 


12-28 


14-27 


27-33 


8-32 


12-28 


10-30 


26-33 


Al 


-1.02(17) 


-1.09(20) 


-1.01(12) 


-0.49(15) 


-0.92(12) 


-0.97(19) 


-0.89(10) 


-0.38(14) 




0.3 


4.3 


2.9 


3.0 


1.7 


4.2 


1.8 


2.1 




12-28 


10-30 


12-28 


28-33 


10-28 


13-27 


14-27 


25-34 


Al 


0.01(5) 


-0.01(3) 


-0.03(3) 


0.12(7) 


-0.01(4) 


0.02(5) 


-0.01(3) 


0.14(9) 




3.9 


4.9 


0.9 


1.8 


1.2 


4.0 


1.1 


0.9 




12-28 


12-30 


11-29 


26-33 


12-24 


14-28 


11-29 


26-33 


T/ 


5.66(69) 


5.09(59) 


5.27(30) 


2.73(52) 


4.90(91) 


4.79(66) 


4.75(39) 


2.76(47) 




4.0 


2.9 


0.7 


2.3 


2.5 


3.1 


1.3 


1.5 




12-28 


13-30 


12-28 


26-33 


13-25 


14-28 


10-30 


26-33 


%' 


5.86(77) 


5.39(49) 


5.38(32) 


3.24(48) 


4.87(121) 


5.01(46) 


4.93(30) 


2.96(48) 



Table 3a. The same as in Table 2a, but for k = 0.155. 
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PS{^\SP 


PW{^\WP 


PS^\WP 


PS\\\WP 


AS{^\SA 


Aiy[o]iyA 


AS^\WA 


AS{\\WA 




2.6 


2.7 


0.6 


1.0 


2.0 


1.6 


0.9 


0.8 




10-30 


13-29 


10-30 


24-32 


12-28 


14-28 


11-29 


26-33 


^2 


-19.7(26) 


-18.4(15) 


-18.0(11) 


-9.0(15) 


-16.7(23) 


-16.8(13) 


-15.8(12) 


-9.0(11) 




2.7 


4.0 


2.4 


1.3 


2.2 


3.3 


1.1 


2.2 




10-30 


13-27 


10-30 


26-35 


12-28 


14-28 


15-27 


26-33 


5^ 


-1.45(32) 


-1.44(33) 


-1.43(17) 


-0.94(23) 


-0.98(31) 


-1.24(14) 


-1.05(24) 


-0.82(15) 




5.2 


3.9 


1.0 


1.4 


1.6 


2.5 


1.0 


1.6 




11-29 


14-28 


12-26 


26-33 


14-27 


16-30 


12-28 


27-32 


v.^ 


0.07(4) 


0.06(4) 


0.08(2) 


0.00(4) 


0.13(11) 


0.02(7) 


0.05(4) 


-0.08(9) 




3.8 


3.4 


2.0 


1.4 


4.6 


4.3 


1.7 


1.3 




12-28 


14-27 


10-30 


26-33 


14-27 


14-28 


15-27 


26-31 


v^ 


0.04(1) 


0.04(1) 


0.04(1) 


0.01(1) 


0.02(3) 


0.04(2) 


0.03(1) 


0.03(2) 




7.4 


3.6 


2.5 


1.8 


2.8 


1.5 


1.4 


1.3 




12-28 


11-31 


10-30 


26-35 


14-27 


15-27 


11-29 


26-34 


Al 


-1.06(17) 


-0.96(6) 


-0.97(9) 


-0.52(9) 


-0.89(14) 


-0.86(8) 


-0.87(13) 


-0.42(8) 




0.1 


3.9 


1.9 


1.2 


0.3 


5.9 


1.4 


1.1 




12-28 


15-25 


16-24 


26-33 


12-28 


14-28 


12-25 


28-34 


Ai 


0.57(10) 


0.51(9) 


0.49(9) 


0.61(13) 


0.49(10) 


0.48(11) 


0.45(4) 


0.53(8) 




4.0 


3.3 


2.6 


2.2 


4.2 


3.8 


4.1 


2.1 




15-25 


13-27 


14-27 


28-34 


12-28 


14-28 


9-31 


26-33 


Tl 


0.14(6) 


0.13(3) 


0.14(3) 


0.03(2) 


0.09(3) 


0.11(3) 


0.11(4) 


0.06(3) 




4.4 


1.6 


2.6 


1.0 


1.9 


3.0 


1.8 


0.8 




14-27 


14-27 


12-28 


26-33 


14-26 


17-26 


10-30 


26-35 


T? 


0.15(4) 


0.14(1) 


0.14(2) 


0.05(2) 


0.06(5) 


0.12(8) 


0.13(4) 


0.07(2) 



Table 3b. The same as Table 3a, but for the two color loop contribution. 
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K = 0.154 


K = 0.155 


alff 


lia 


Bk(P = 0) 
12-28 


BkiP = l) 
24-34 


Bk 


Bk(P = 0) 
10-30 


BkiP = l) 
26-34 


Bk 


0.00 


— 


0.298(44) 


0.446(42) 


0.83(21) 


- 0.023(63) 


0.263(62) 


0.77(19) 


1.00 


1.0 


0.373(45) 


0.494(52) 


0.81(23) 


0.109(74) 


0.330(61) 


0.71(20) 


1.00 


TT 


0.362(42) 


0.476(50) 


0.78(22) 


0.110(71) 


0.319(60) 


0.69(22) 


1.00 


1.77r 


0.356(42) 


0.468(50) 


0.76(22) 


0.110(70) 


0.315(59) 


0.67(21) 


1.338 


1.0 


0.405(44) 


0.512(55) 


0.79(24) 


0.175(75) 


0.359(69) 


0.68(24) 


1.338 


TT 


0.388(42) 


0.486(53) 


0.75(22) 


0.176(71) 


0.344(62) 


0.64(22) 


1.338 


1.77r 


0.380(41) 


0.475(51) 


0.72(22) 


0.176(69) 


0.337(61) 


0.62(22) 


1.75 


1.0 


0.453(45) 


0.536(58) 


0.75(23) 


0.268(70) 


0.401(74) 


0.64(25) 


1.75 


TT 


0.429(42) 


0.497(53) 


0.68(22) 


0.268(63) 


0.379(71) 


0.57(23) 


1.75 


1.77r 


0.417(41) 


0.497(51) 


0.64(21) 


0.268(59) 


0.369(69) 


0.55(22) 



Table 4. The lattice S— parameter for the perturbatively improved operator O for zero 
and one unit of lattice momentum, and the S-parameter after momentum subtraction. 
The data show the magnitude of the variation with the value of g'^g and //a used in the 
perturbative renormalization constants. The range of time-slices used in the fits is specified 
in the header. 
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Wilson 
K = 0.154 


Staggered 
niq = 0.02 + 0.03 


Wilson 
K = 0.155 


Staggered 
niq = 0.01 + 0.02 


Mk 


0.370(6) 


0.374(3) 


0.297(11) 


0.296(2) 


Ek 


0.511(12) 




0.466(22) 




VI 


-0.57(33) 


-0.48(2) 


-0.71(53) 


-0.79(4) 


v^ 


-0.11(8) 


-0.036(4) 


-0.13(11) 


-0.10(1) 


VI 


-0.10(7) 


-0.056(3) 


-0.12(15) 


-0.13(1) 


VI 


0.01(2) 


-0.009(1) 


-0.02(4) 


-0.024(2) 


A\ 


0.28(18) 


0.15(1) 


0.41(42) 


0.36(2) 


Al 


0.27(16) 


0.063(5) 


0.29(32) 


0.13(1) 


A\ 


0.28(12) 


0.33(1) 


0.40(20) 


0.42(1) 


A\ 


0.79(26) 


0.81(1) 


0.82(40) 


0.85(2) 


Ba 


1.62(42) 


1.35(3) 


1.92(89) 


1.76(5) 


Bv 


-0.78(40) 


-0.58(2) 


-0.98(71) 


-1.04(5) 


Bk 


0.83(21) 


0.76(1) 


0.77(19) 


0.72(2) 



Table 5. Comparison of individual S— parameters, for space/time and l-loop/2-loop 
components of the vector and axial four-fermion operators, between Wilson and staggered 
fermions at matching values of the kaon mass. The kaon energy at p = (0, 0, 0) and 
p = (0, 0, 1) measured from the 2-point correlators is also given. All results are quoted for 
9!s = 0. 
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5i(o) 


BkW 


Mk 


E{1) 


Bk 


glff = 0.000 

K = 0.159 g^ff = 1.091 
glff = 1-909 


0.327(88) 
0.389(74) 
0.401(56) 


0.507(96) 
0.537(80) 
0.505(58) 


0.477(13) 


0.660(16) 


0.98(43) 
0.92(36) 

0.78(27) 


glff = 0.000 

K = 0.160 glff = 1.091 

glff = 1-909 


-0.593(88) 
-0.415(76) 
-0.217(57) 


-0.052(51) 
0.051(44) 
0.132(34) 


0.362(7) 


0.562(22) 


0.93(24) 
0.89(21) 
0.76(17) 


<72^.^ = 0.000 

K = 0.156 glff = 1.071 

^e'// = 1-875 


0.497(33) 
0.534(31) 
0.512(26) 


0.685(39) 
0.693(36) 
0.629(30) 


0.456(5) 


0.613(10) 


1.23(20) 
1.15(19) 
0.97(16) 


glff = 0.000 

K = 0.157 glff = 1.071 

glff = 1-875 


0.154(25) 
0.161(25) 
0.234(21) 


0.436(29) 
0.438(29) 
0.439(24) 


0.358(7) 


0.501(16) 


1.14(19) 
1.13(19) 
0.95(16) 



Table 6. The same as in Table 4, but for lattices generated with two flavors of dynamical 
Wilson quarks. The upper and lower halves of the table correspond to /3 = 5.5 and /3 = 5.6 
respectively. These numbers are obtained with //a = 1.0. 
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K = 0.154 


K = 0.155 


alff 


(la 


^3/2 


or 


^3/2 


03/2 


0.00 


— 


0.902(24) 


0.947(30) 


0.908(44) 


0.963(50) 


1.00 


1.0 


0.871(24) 


0.918(29) 


0.878(44) 


0.935(50) 


1.00 


TT 


0.918(29) 


0.945(30) 


0.911(46) 


0.962(52) 


1.00 


1.77r 


0.901(25) 


0.953(30) 


0.921(46) 


0.969(52) 


1.338 


1.0 


0.832(23) 


0.877(23) 


0.830(43) 


0.894(48) 


1.338 


TT 


0.891(25) 


0.934(30) 


0.902(46) 


0.951(51) 


1.338 


1.77r 


0.908(25) 


0.950(30) 


0.921(47) 


0.966(52) 


1.75 


1.0 


0.747(21) 


0.785(26) 


0.753(40) 


0.802(45) 


1.75 


TT 


0.869(25) 


0.911(30) 


0.883(46) 


0.928(51) 


1.75 


1.77r 


0.901(26) 


0.941(31) 


0.916(48) 


0.958(52) 



Table 7. The S-parameter for the LR electromagnetic penguin operators on the quenched 
lattices. The fit range is t = 13 — 27 in all cases. The results are shown for a number of 
values of gl^^ and (la used in the perturbative renormalization constants. 
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